Breaking chirality in nonequilibrium systems on the lattice 
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Abstract. - We study the dynamics of fronts in parametrically forced oscillating lattices. Using 
as a prototypical example the discrete Ginzburg-Landau equation, we show that much information 
about front bifurcations can be extracted by projecting onto a cylindrical phase space. Starting 
from a normal form that describes the nonequilibrium Ising-Bloch bifurcation in the continuum 
and using symmetry arguments, we derive a simple dynamical system that captures the dynamics 
of fronts in the lattice. We can expect our approach to be extended to other pattern-forming 
problems on lattices. 
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Extended systems on lattices have played a major role 
in the development of nonlinear science. One may recall 
celebrated models as the discrete sine-Gordon equation or 
the Fermi-Pasta-Ulam experiment. Usually the develop- 
ment of the field has been associated with conservative 
systems, but dissipative lattices have attracted growing 
attention in the last two decades (see [1] for a review). 
Two prominent examples of such lattices are provided by 
the discrete version of the Nagumo and Ginzburg-Landau 
partial differential equations. The former has been pro- 
posed as a model of myelination of neuronal fibres [2] ; and 
the latter describes, among others, dissipative solitons [3], 
the dynamics of lines of vortices [4] and coupled wakes [5] 
in hydrodynamics. 

The complex Ginzburg-Landau equation [6] universally 
describes the dynamics of an extended medium in the 
neighbourhood of an oscillatory instability. Under homo- 
geneous resonant n : 1 forcing, different regions of space 
may lock to the driving with different phase relations; and 
domain walls separate these regions. The prototypical 
2 : 1 resonant case leads to the so called parametrically 
forced complex Ginzburg-Landau (FCGL) equation. This 
equation has been subject of extensive study since the sem- 
inal work by GouUet and coworkers in ref. [7]. There it 
was found a front bifurcation which is the nonequilibrium 
analogue of the Ising-Bloch transition in ferromagnets. 

In this Letter, we show that the dynamics of fronts in 
the FGGL equation on the lattice is captured by a nor- 
mal form consisting of two ordinary differential equations. 



The bifurcations linking different dynamics of the front 
(including bistable regimes) are observed both in a pro- 
jection of the system's variables onto a cylindrical phase 
space and in the normal form. Our results are relevant for 
experiments where discretisation is given as in arrays of 
coupled pendula [8,9], electronic circuits [10,11], or chem- 
ical systems [12]; and also in systems usually modeled as 
continuous, but that are intrinsically discrete (or behave 
like a lattice due to a spatially periodic modulation of the 
medium) . 

For a lattice, the FCGL equation [13,14] takes the form: 

Aj = (1 + tiy)Aj - (1 + if3)\Aj\^Aj + jA* 

+K{l + ia){A.j+i+Aj-i~2Aj), (1) 

where Aj = pje^^^ is a complex variable. The parameter 
7 measures the forcing strength, and k controls the cou- 
pling between neighboring units. Parameters v, (3, and 
a account for the detuning, the nonisochronicity, and the 
dispersion, respectively. 

The continuum limit. — First of all we recall the 
results for the continuous version of (1): 

dtA ^ {l + iv)A- {\ + il3)\A\'^ A + -fA* + {l + ia)l^A . (2) 

Vanishing values of v, f3, and a allow to cast (2) into a 
variational form: dA/dt = —5T/5A*. In this case, stable 
front solutions minimise the energy functional J^, and de- 
pending on the forcing 7 they can be chiral (7 < 7/_b = 
1/3) or achiral (7 > 7/5). If A vanishes at the centre of 
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Fig. 1: Regions in the parameter space iz-'y where Ising (I) and travelling or stationary Bloch fronts (denoted by TB and SB, 
respectively) appear: (a) in the continuum, and (c) on the lattice, k = 0.25 in eq. (1). Sketches of all possible fronts (b) in the 
continuum, and (e) on the lattice. Whereas the fronts I, TB and SB in (e) are stable for some parameter values, the front U is 
always unstable. 



the front the so-called Ising front is found, otherwise the 
front is chiral and A does not vanish anywhere: two such 
Bloch solutions with opposite chirality exist. 

The nonequilibriuni Ising-Bloch (NIB) transition is ob- 
served for nonzero values of i^, /3, or a [7,15]. The terms v, 
P and a are perturbations to a gradient system, and cause 
the Bloch fronts to move, in either direction depending on 
their chirality (see [14] and references therein). 

Front dynamics. — We restrict ourselves hereafter 
to the case (3 — a — 0, so that v remains as the pa- 
rameter that breaks the variational character of the sys- 
tem; we note nevertheless that the same qualitative results 
are obtained perturbing variationality through ^, /3, and 
a. Inside the region 7 > \v\ (the Arnold tongue), see 
fig. 1(a), the local dynamics is bistable with two stable 
fixed points As± = ipse*'''^ with ps = [l + Vl^ - t^^V^^, 
and ipS G (— 7r/2,7r/2) the solution of sin(2?/'s) — i^/l, 
cos(2-05') — {pg — The fixed point at the origin 

= (0,0) is either completely unstable (7^ < 1 + v^), 
or saddle (7^ > 1 + ly^). We are interested in the dy- 
namics of fronts connecting the two stable fixed points: 

In the contimmm, the locus of the NIB transition can 
be calculated analytically [16]: "fNiB{v) = [vT+9j^]/3, 
see fig. 1(a) [and fig. 1(b) for a sketch of both front types]. 
As mentioned above, for v — Q we recover the variational 
case, ^nib{i^ = 0) — ^ib, and Bloch fronts are stationary. 

Our extensive numerical simulations of the FCGL equa- 
tion on the lattice^ have revealed several front types not 
present in the continuum. These fronts are sketched in 

^The numerical calculations were performed with zero flow 
boundary conditions Aq = Ai, j4]v+i = Ajv, and a lattice size 
N large enough to neglect boundary efi'ects on the front dynamics 
(with N typically being 128). 



fig. 1(e) and figs. l(c,d) show their regions of stabil- 
ity. Stationary Bloch fronts (SB) exist on a finite region 
of parameter space around an interval of the line 1^ = 0. 
Additionally, two regions of bistability are found: In one 
of them Ising and travelling Bloch fronts coexist (I-I-TB). 
The other bistable region (SB+TB) is shown in fig. 1(d), 
and it is a small triangle with stable stationary and trav- 
elling Bloch fronts. 

Cylindrical coordinates. — A projection of the 2N 
degrees of freedom onto a two-dimensional phase space 
greatly simplifies the analysis of the fronts. One way of 
performing this projection is 

N 

$ = Re[y ^e'('^^-'^^)], (3) 

N 

C = im[y -^e*(^-'^s)] . (4) 

This corresponds to a rotation and compression of the 
complex plane A such that the stable fixed points are now 
located on the real axis at ±1. This projection permits 
us to discern if a front is symmetric with respect to the 
origin (in such a case C = 0). Note that $ is a cyclic vari- 
able that takes the same value when the front advances 
or recedes one lattice unit (provided the front is far from 
the boundaries) and thus can be defined modulo 2. Vari- 
able C measures the deviation from stationarity and is 
intrinsically bounded. 

Typical examples of the dynamics in the reduced coor- 
dinates (3)-(4) are shown in figs. 2(a,b). Two stationary 
fronts, both with C = 0, exist for all parameter values. 
One, the Ising front (I), is located at C) = (0, 0) if the 
number of units N is even (or at (1,0) if N is odd) and 
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Fig. 2: Numerically obtained flows projecting the front dynam- 
ics onto the cylindrical phase space $-C defined by eqs. (3) and 
(4). In (a) an example of coexistence I+TB is shown (i^ = 0.1 
and 7 = 0.183). In (b) SB+TB {f = 0.045 and 7 = 0.148). 
Symbols •, x , o denote stable, saddle and unstable fixed points, 
respectively. 

is a continuation for k > of the trivial solution at zero 
coupling: [. . . , As+, As+, As+, As_, As-, As-, . . .]. The 
second stationary solution, denoted U, is at (1,0) (respec- 
tively, at (0, 0) for odd N), and is a continuation of the un- 
stable front solution: [. . . , ^5+ , ^5+ , Aq , As- ,As-,.-.]- 
For some parameter values, there exist extra fixed points 
located off the # axis. We label these chiral solutions sta- 
tionary Bloch (SB) fronts. Travelling Bloch (TB) fronts 
correspond to periodic orbits around the cylinder. Due 
to symmetry they always appear in pairs circulating in 
opposite directions. 

Normal form. — Next, wc present a simple ordinary- 
differential-equation model that generates dynamics like 
the front dynamics on the FCGL lattice. We will obtain 
this normal form via symmetry arguments. 

In continuous systems the NIB transition is a parity- 
breaking (pitchfork) bifurcation coupled to a translation- 
invariant coordinate. It is an example of the so-called 
drift-pitchfork bifurcation found in a number of situations 
(e.g. [17]), usually as a secondary instability. Its normal 
form is [18, 19] (f> = c, c = {jj, — c^)c, where /i is the bifur- 
cation parameter, e.g., yit oc {jnib — 7). Coordinates </>, 
c represent the position and the velocity of the concerned 
structure (the front in our case). 

Knobloch et al. [20] considered the breakdown of the 
continuous translational invariance to study the parity 
breaking of a periodic pattern in the presence of an in- 
homogeneity. This leads to the inclusion of small periodic 
terms (sinusoidals in the simplest case) that preserve the 
invariance under inversion (0, c) — > {—^, — c): 

^ = c—esm<p, (5) 
c = (/i + (5cos0 — c^)c + 7? sin 0, (6) 

where (f) is now an angular variable. In [21] it was sug- 
gested that this normal form could also be used to anal- 
yse parity-breaking bifurcations found in discrete bistable 



media {e.g. arrays of FitzHugh-Nagumo units). 

In the continuum, as the variational limit is approached 
the velocity of the Bloch front decreases (and becomes zero 
at variationality). Consequently, we introduce a small pa- 
rameter X that accounts for the deviation from the vari- 
ational case. For the situation considered here, x is pro- 
portional to^ v. We have then at leading order: 

<P = xc-esm<j>, (7) 
c = {n + S cos (j) — c'^)c . (8) 

The term rj sin (p present in (6) is absent in (8) due to the 
invariance of the discrete FCGL equation under the trans- 
formation {ly, f3, a. A) {—iy,—l3,—a,A*). This requires 
the symmetry under (x, (j), c) (— x, 0, — c) to be satisfied. 
This considerably simplifies the normal form. 

Figure 3(a) shows the regions of the parameter space 
X-M) and fig. 3(c) is a sketch of the corresponding phase 
spaces. The normal form (7)-(8) qualitatively reproduces 
front dynamics, which can indeed be projected on a cylin- 
der, e.g. the phase space $-C introduced above. There 
are two pitchfork bifurcations a.t n = ±S and the pa- 
rameter space is organised by two codimcnsion-2 points. 
We may see in fig. 3(b) that at a degenerate pitchfork 
point located at x^ = — 2e^/(5 the pitchfork bifurcation 
line switches between subcritical and supercritical. The 
second codimension-2 point is a saddle-node separatrix- 
loop point [22] where the saddle-node bifurcation on the 
invariant circle (SNIC) splits into an off-cycle saddle-node 
{SN) and a homoclinic {Horn) bifurcations. 

In contrast to the continuous case, travelling Bloch 
fronts appear with nonzero chirality, and the (average) 
velocity of the front v asymptotically follows: (i) the 
familiar square-root law when crossing the SNIC line; 
(ii) a logarithmic law below the homoclinic line: = 
a — {l/\u) ln(7 — jHom), where A„ is the positive eigen- 
value corresponding to the saddle Bloch front. 

Let us finally discuss the parameters modelling the dis- 
cretisation strength: e and 5. In the continuum limit, our 
numerical simulations indicate that discretisation enters in 
the normal form through S at order 0{k~^), whereas e van- 
ishes much faster (possibly exponentially) as k ^ 00. In 
addition, assuming opposite signs for e and S correctly in- 
herits the bifurcations in the anticontinuum limit (k ^ 0). 
Accordingly, our choice e = 0.1 and S = —0.2 in fig. 3(a) 
satisfies these constraints {i.e. \e\ < \5\, eS < 0). The or- 
ganisation of the parameter space is qualitatively robust 
to changes of both parameters in a wide range. 

Phaise equation. — Normal form (7)- (8) is obtained 
in a perturbative way including small terms that model 
breakdown of translational invariance, and the departure 
from {v,(3,a;^) = {0;^ib), the equilibrium Ising-Bloch 
bifurcation point. This means that we cannot predict the 
behaviour of the bifurcation lines far from this point. In 

^According to Eq. (4) in [7], in the Vciriational limit one has x oc 

u-(3-{a-P) -t. 
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Fig. 3: (a,b) Parameter space of the normal form (7)-(8) with 
e = 0.1 and 5 = —0.2; the loci of local and global bifurcations 
are depicted with solid lines. The black dots in the inset (b) 
indicate the location of two codimension-2 bifurcation points, 
(c) Schematic representation of the phase space (p-c for the 
five regions in (a,b). Gray-dashed lines indicate the link be- 
tween different states via codimension-1 bifurcations {PF and 
PFsub are, respectively, super- and subcritical pitchfork, and 
SN, Horn, and SNIC stand for (off-cycle) saddle-node, homo- 
clinic, and saddle-node on an invariant circle bifurcations). 



particular, the (SNIC) bifurcation line, which separates 
SB and TB regions in fig. 1(c), can be proved to end at 
= 7 = as results from the following facts. For small 
7, a phase reduction [7, 14] of the FCGL lattice yields a 
discrete Nagumo-type equation [the overdamped Frenkel- 
Kontorova (FK) model]: = v — ^sin[2'ipj) + «;(f/'j-i-i + 
il^j-i — '^ipj)- Interestingly, v acts as a symmetry-breaking 
parameter, but the symmetry of the original model is hid- 
den in the factor 2 inside the sm{2ipj) term, which allows 
two mirror front solutions connecting ipj^-oo = V'S a-nd 
ij^j^oo = V's + Discreteness implies the existence of an 
interval of 'propagation failure' [2] where propagation is 
blocked for nonvanishing v. Well established results for 



the FK model (see e.g. [23,24]) state that the propagation 
threshold {SNIC line) should approach the 7 axis expo- 
nentially fast^ as 7 ^ 0. 

Conclusions. — The nonequilibrium Ising-Bloch bi- 
furcation in the parametrically driven complex Ginzburg- 
Landau equation is one of the most studied pattern in- 
stabilities. In the continuum, breaking of chirality causes 
the front to move. However, as shown here, on the lattice 
a more complex scenario appears: specifically, two types 
of bistability and a region with chiral (Bloch) stationary 
fronts. We have demonstrated that a normal form with 
two variables captures the dynamics of the front and the 
bifurcations between different regions. It is to be empha- 
sised that the normal form (7)-(8) is based on symmetry 
arguments that provide general qualitative results inde- 
pendent of details as, for instance, the parameter of non- 
variationality or the discretisation-order of the Laplacian. 

The dynamics of patterns on lattices is typically much 
harder to solve analytically than in their continuous coun- 
terparts. The continuum, usually serves as a zeroth order 
approximation that is not necessarily accurate. Discrete- 
ness typically introduces new dynamical regimes as shown 
in the current letter for one spatial dimension*. Through 
a modification of the normal form for the continuum, we 
have been able to reproduce the dynamics of fronts on 
a discrete medium and the structure of parameter space. 
This approach should also work in other problems on lat- 
tices. 
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